Rayleigh scattering, mode coupling, and optical loss in silicon microdisks 



O 

o 

(N 



O 

• I— I 
-(— > 

Oh 

q 
o 

a.: 



> 



\0' 

o- 

^: 

o, 

o : 

•1— I 

Oh. 



X 



Matthew Borselli, Kartik Srinivasan, Paul E. Barclay, and Oskar Paintecl 
Department of Applied Physics, California Institute of Technology, Pasadena, CA 91125, USA. 

(Dated: February 2, 2008) 

High refractive index contrast optical microdisk resonators fabricated from silicon-on-insulator 
wafers are studied using an external silica fiber taper waveguide as a wafer-scale optical probe. 
Measurements performed in the 1500 nm wavelength band show that these silicon microdisks can 
support whispering-gallery modes with quality factors as high as 5.2 x 10^, limited by Rayleigh 
scattering from fabrication induced surface roughness. Microdisks with radii as small as 2.5 /im are 
studied, with measured quality factors as high as 4.7 x 10^ for an optical mode volume of 5.3(A/n)^. 

PACS numbers: 42.60.Da, 42.81. Qb 



Recent studies of optical resonators in glass based 
microspheres m 12, |3] I rnicrorings 4] , and microtoroids 5| 
have highlighted the applications afforded by the ex- 
tremely long photon lifetime of whispering-gallery-modes 
(WGMs) supported by these structures. Furthermore, 
recent work by Ilchenko et al. illustrated the advantages 
of creating WGMs in crystalline materials[S|. Applica- 
tions for such devices include quantum networking 7], 
low threshold non-linear optical sourcesjSJ, and compact 
micro-optical circuits[j|. The ability to create similar 
high quality factor (Q) WGM resonators in III-V or sil- 
icon (Si) semiconductors has thus far been hampered by 
the large refractive index of most semiconductors and the 
resulting sensitivity to surface roughness |9|, [lOJ. In this 
Letter we describe measurements of micron-sized Si mi- 
crodisk resonators supporting transverse-magnetic (TM) 
WGMs with significantly reduced sensitivity to disk-edge 
roughness. These modes have measured Q values as high 
as 5.2 X 10^ and effective modal volumes (T4ff) as small as 
5.3 cubic wavelengths in the material. The largest Q/Vcs 
ratio is measured to be 8.8 x 10^, greater than the val- 
ues measured in ultra-small volume photonic crystals [llj 
and comparable to the values measured in ultra-high-Q 
microspheres and microtoroids ^^^]. 

The silicon microdisks in this work are fabricated 
from a silicon-on-insulator (SOI) wafer consisting of a 
344 nm thick p-doped Si layer of resistivity 1-3 rj-cm 
atop a two micron Si02 layer. Processing of the mi- 
crodisks begins with the deposition of a 20 nm Si02 
protective cap layer using plasma-enhanced chemical- 
vapor-deposition. Electron beam lithography is used 
to create a polymer resist etch mask, and a low- 
bias voltage inductively-coupled-plasma reactive-ion-etch 
with SF6:C4Fg gas chemistry [llj then transfers the cir- 
cular microdisk pattern into the top Si layer. After dry- 
etching, the sample is immersed in buffered hydrofluoric 
acid to undercut the bottom Si02 cladding, as shown in 
Figure n The thin 20 nm Si02 top cap layer is also re- 
moved in this process, providing a clean, smooth top Si 
surface. A final rinse in deionized water is performed, 
followed by a high-purity nitrogen spray drying step. 

In order to characterize the microdisk resonators, 



an evanescent fiber taper coupling technique is 
employed[ll|. In this process, an optical fiber is adiabat- 
ically drawn to a 1-2 /j,m diameter so that its evanescent 
field is made accessible to the environment. In this work, 
the fiber taper is positioned to the side of the microdisks, 
with a center height equal to that of the middle of the 
microdisk. Measurements of the taper transmission as a 
function of the lateral taper-microdisk gap {g) are then 
performed using a swept wavelength tunable laser source 
(A = 1509-1625 nm) with fine frequency resolution of 10 
MHz. A set of paddle wheels are used to adjust the po- 
larization state of the fiber taper mode in the microdisk 
coupling region, providing selective coupling to the TE- 
like (TM-like) WGMs with dominant electric field paral- 
lel (normal) to the plane of the microdisk. For the 344 
nm Si layer thickness of the microdisks studied here, only 
the fundamental vertical slab mode for TE and TM po- 
larization are strongly guided. As such, we only consider 
the fundamental vertical modes of the microdisks in what 
follows, labeling them simply by Pm.m where p is either 
TE or TM, and m and n are the characteristic radial and 
azimuthal number, respectively [l^- 

Microdisks of two different sizes, radius i? = 2.5 and 
4.5 /J,m, are fabricated and tested. A broad wavelength 
scan covering the 1509-1625 nm wavelength range is 
initially employed to map out the different microdisk 
modes. As wiU be detailed elsewhere [l3|. the adjustable 
polarization state in the taper along with the WGMs' 
strength of coupling and linewidth is used to determine 
sets of modes with a common free spectral range. These 
measurements provided a reliable determination of ra- 
dial mode number. An effective index two-dimensional 
model is then used to estimate the azimuthal number. 
Using this mode identification technique we found that 
the highest Q modes in both sizes of microdisks are con- 
sistently of TM polarization, and corresponded to the 
lowest radial number, n ^ 1. 

The inset of Fig. ^a.) shows the evanescent coupling to 
a TM44_i WGM of a i? = 4.5 /im microdisk, with tapered 
fiber positioned 1.1 /im laterally from the disk edge. The 
observed double resonance dip (doublet) is a result of 
Rayleigh scattering from disk surface roughness as dis- 




-2 2 

arclength (\im) 



^6 

CD 



M A 


— data 
""gaussian fit 


/ \ 


- 





-600 



-400 



-200 200 

displacement (nm) 



400 



600 



FIG. 1: SEM micrographs of a _R = 2.5 ^m Si microdisk: (a) 
side-view illustrating Si02 undercut and remaining pedestal, 
(b) high contrast top-view of disk, and (c) zoomed-in view 
of top edge showing disk-edge roughness and extracted con- 
tour (solid white line), (d) Plot of extracted contour versus 
arclength. (e) Autocorrelation function of the microdisk con- 
tour and its Gaussian fit. 



cussed below, which lifts the degeneracy of clockwise 
{cw) and counter-clockwise {ccw) propagating WGMs in 
the microdisk[2|. Fitting the shorter wavelength mode 
of the doublet to a Lorentzian yields a loaded linewidth 
of 3.9 pm with a 5% coupling depth. These measure- 
ments are repeated for varying taper-microdisk gaps and 
are recorded in Figure|3a,b). For g > 0.63 /zm, the data 
follows a two-port coupled mode theory with simple ex- 
ponential loading dependence on taper-microdisk gap[2|. 
Fits based upon this model are shown as a solid line in 
each of the plots of Fig. |21 The fiber loading of the mi- 
crodisk is characterized here by a dimensionless effective 
quality factor, Qfibor (inset to Fig. IHb)). The asymptotic 
unloaded linewidth is found to be 3.0 pm for this WGM, 
corresponding to a bare-cavity Q of 5.2 x 10^. The ob- 
served non-exponential loading dependence for small g is 
due to the large phase mismatch between the glass fiber 
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FIG. 2: Fiber taper measurements of a TM44,i WGM of a 
microdisk with R = 4.5 fim. (a) Lorentzian full-width half- 
maximum (FWHM) linewidth versus taper-microdisk gap. 
(inset) Taper transmission showing high-Q doublet, (b) Res- 
onant transmission depth versus taper-microdisk gap. (inset) 
loading versus taper-microdisk gap. 
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(2.5,TM23,i) 


5.3 
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12.2 


(4.5,TE5o,i) 


16.2 



4.7 (3.5) X lO'^ 93 (112) 

9.0 (9.6) X 10"* 126 (146) 

5.2 (6.2) X 10^ 37 (46) 

1.7 (1.8) X 10^ 78 (59) 



TABLE I: Summary of theoretical and measured mode pa- 
rameters for R = 2.5 and 4.5 fim Si microdisks. Theoreti- 
cal surface-scattering values for Qss and A\±m are shown in 
parentheses. 



and Si microdisk modes, further studies of which will be 
presented elsewhere jl3l|. 

Similar measurements are performed for all n ~ 1 
modes (of both polarizations) in each of the two different 
microdisk sizes, and a summary of the measured bare- 
cavity Q and doublet mode-splitting values are given in 
Table m In order to understand the limiting loss mech- 
anisms and the observed trend of Q with polarization 
and radial mode number, a complementary theoretical 



analysis is also performed. Both the radiation limited 
Q (^ 10^^[lJ|) and the free carrier absorption Q of the 
pooped Si disk layer (estimated to be > 10'^ from Ref. 
[15|), are far greater than the measured Q values. Sur- 
face state absorption is another possible loss mechanism; 
however, several trials are performed with differing fi- 
nal chemical cleaning steps with no perceivable change 
in Q. The rather large measured splitting of the cw 
and ccw traveling wave modes indicates strong surface 
scattering [la|, and therefore we believe the dominant 
source of optical loss to be scattering from index per- 
turbations at the microdisk surface. 

Surface roughness on the top and bottom surfaces of 
the microdisks analyzed in this work is negligible in com- 
parison to the azimuthal variation of the disk radius 
[Fig. n^c)], a result of mask erosion and deformation 
during the dry etching step. For this type of rough- 
ness, the index perturbation can be approximated by 
6e = eoSn'^hAr{s)S{r — R)5{z), where Eq is the free space 
permittivity, Sn^ — n'^ — n'^, Ud ^ 3.46 is the Si refractive 
index, rio ~ 1 is the index of the surrounding air, h is the 
disk height, and Ar(s) is the radial surface roughness rel- 
ative to the unperturbed disk radius,17j. The parameters 
r, s, and z correspond to the radius from the disk center, 
arc length along the disk perimeter, and height along the 
disk edge, respectively. These perturbations of the disk 
radius set up polarization currents, J = —iuJoSeE, which 
couple the WGMs of the perfect microdisk to radiation 
modes or other nearly resonant WGMs. The Volume 
Current Method[l^ can be used to calculate the radi- 
ated power into free space from J, providing an estimate 
for the surface-scattering quality factor, Qss- In the limit 
that the correlation length of the roughness (Lc) is much 
smaller than the wavelength in the material. 
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where V = ^/RL^harSn^ , Us{e) is the average e-polarized 
electric field energy density at the disk edge normalized 
to unit modal energy, ar is the standard deviation of 
the roughness, and G(e) = {4/3, 2, 2/3} is a geometri- 
cal radiation factor for the e = {r, 0, z} different electric 
field polarizations. The mode coupling between counter- 
propagating cw and ccw traveling wave modes can also 
be found via a time-dependent perturbation theory [ig. 
This coupling breaks the degeneracy of the cw and ccw 
modes (±m), creating two standing wave modes with 
wavelength splitting 
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which edge contours with near nanometer resolution can 
be obtained [Fig. H^d)]. From a Gaussian fit to the auto- 
correlation function of the surface roughness [Fig. H^e)], 
estimates for ar and L^, are found to be 2.8 nm and 40 nm, 
respectively, typical of the fabricated microdisks in this 
work. Using the results from this roughness analysis in 
Eqs. QJ and l^l, theoretical estimates for Q and AA±m 
are calculated and tabulated in Tabled with good agree- 
ment found between theory and the measured microdisk 
properties. In these calculations, an effective index model 
based on Ref. |l3| is used to estimate Us{e) and the ef- 
fective mode volume Voff|i2ji with the approximate TM 
modes given by a single z field component and the TE 
modes consisting of f and (f> electric field components. 
The salient result of this analysis is that a TE mode has 
roughly three times the radiated power versus a compa- 
rable TM mode. This is explained by the fact that both 
the f and components couple to radiation modes more 
strongly than the z component of the field. The modes 
of higher n number are also measured and found to be 
lower Q as predicted by the disk-edge surface-scattering 
theory, although for n ^ 2 the effects of the underlying 
Si02 pedestal may be significant. 

In summary. Si microdisks of a few microns in radius 
have been fabricated with measured Q > 5 x 10^. In con- 
trast to theoretical predictions for an ideal microdisk, we 
find that the TM modes have significantly higher mea- 
sured Q values than TE modes due to their inherent 
reduced sensitivity to disk-edge surface roughness. In 
comparison to previously measured high-Q semiconduc- 
tor microdisks |^ , the Q of the smallest WGMs measured 
in this work are > 20 times larger while maintaining a 
similar effective mode volume. 
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where Us = '^^Us{e)- SEM micrographs such as those 
shown in Figs. n^b,c) provide high contrast images from 
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